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ABSTRACT 

Direct imaging searches have revealed many very low-mass objects, including a 
small number of planetary mass objects, as wide-orbit companions to young stars. 
The formation mechanism of these objects remains uncertain. In this paper we present 
the predictions of the disc fragmentation model regarding the properties of the discs 
around such low-mass objects. We find that the discs around objects that have formed 
by fragmentation in discs hosted by Sun-like stars (referred to as parent discs and 
parent stars) are more massive than expected from the Mdisc ~ relation (which is 
derived for stars with masses M* > 0.2 Mq). Accordingly, the accretion rates onto 
these objects are also higher than expected from the M* — M* relation. Moreover there 
is no significant correlation between the mass of the brown dwarf or planet with the 
mass of its disc nor with the accretion rate from the disc onto it. The discs around 
objects that form by disc fragmentation have larger than expected masses as they 
accrete gas from the disc of their parent star during the first few kyr after they form. 
The amount of gas that they accrete and therefore their mass depend on how they 
move in their parent disc and how they interact with it. Observations of disc masses 
and accretion rates onto very low-mass objects are consistent with the predictions of 
the disc fragmentation model. Future observations (e.g. by ALMA) of disc masses and 
accretion rates onto substellar objects that have even lower masses (young planets 
and young, low-mass brown dwarfs), where the scaling relations predicted by the disc 
fragmentation model diverge significantly from the corresponding relations established 
for higher-mass stars, will test the predictions of this model. 

Key words: Stars: formation, low-mass, brown dwarfs - accretion, accretion discs, 
protoplanetary discs - Methods: Numerical, Hydrodynamics 


1 INTRODUCTION 

Many very low-mass objects, including a small num¬ 
ber of planetary-mass objects, have been observed 
by direct imaging as companions to young stars at 
distances from a few tens to a few hundred AU 


jKrauset_^ Il2008l. I2OI3I: iMa rois et ^ |2008l: |Fahgrtvet_ajJ 


2OO9I: Ireland et al. I2OIII: Kuzuhara et al. 201ll. 

2013; 

Aller et al.l l2013l: iBailev et al.l l2013l: iRameau et al.l 

2013; 

Naud et al.l2014l:lGalicher et aL|2014|). The dominant mech- 


anism for the formation of low-mass stellar and substel¬ 
lar objects (low-mass hydrogen-burning star s, brown dwarfs 
and g iant planets) is st ill uncertain (e.g. IChabrier et al.l 
I2OI4I : IStamatellosI l2014h . It is believed that such ob- 


* E-mail:D.Stamatellos@astro.cf.ac.uk 


jects may form in three ways: (i) by collapsing molec- 
ular cloud cores, i.e. the same way as Sun-like stars 
llPadoan fc Nord lundl |2004 iHennebelle fc Chabrierl l2008l . 
I2OO9I : Ho^kms 2013), (ii) by fragmentation of protostellar 
discs ( |BossI I 1997 I: IStamatelios et al.l 2007al : Attwood et al.l 
I2OO9I : IStamatellos fc WhitworthI 2009al : BolevI 2009ll. which 
may not e ven be centrifugally supported 1 Offner et al.l 
l20ld. l2012ll . and (iii) by ejection o f proto-stellar em¬ 
bryos from their nata l cloud cores llR eipurt h fc Clarl^ 
I2OOII : iBate et"^ l2002l : iGoodwin et alF 2004 1. Addition¬ 
ally, gas giant planets also form by core accretion, i.e. 
by coa gulation of dust particles to progressively larger 
bodies I Safronov fc Zviaginal 19691: Goldreich fc Wa^ll973l: 


Mizund Il980l: Bodenheimer fc PollackI Il986l: Pollack et ahl 


1996ll . Objects formed by core accret ion may even become 


deuterium-burning brown dwarfs fe.g. lMolliere fc Mordasinil 
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| 2012 |) . However, gas giants on wide orbits 100 — 300 AU) 
are not believed to be able to form, at least in-situ, by core 
accretion. 

Substellar objects are difficult to form similarly to Sun¬ 
like stars, and it has been argued that a different mecha¬ 
nism may in fact be at play (e.g. ^SQiitworttL.et_al ] 120071: 
IXhies fc KroupalboOTl : iReggiani fc MeveH 2013l l . A low -mass 
pre-(sub)stellar core has to be very dense and compact 
in order to be gravitationally unstable. Up to now, only 
one clear-cut se h-gravitating brow n dwarf-mass core has 
been observed ll Andre et al.l [ 201 ^ . but such cores have 
small size and they are faint, making them difficult to ob¬ 
serve. Another way to reach the high densities that are 
required for the formation of substellar objects is in the 
discs around young stars. This model has been studied ex¬ 
tensively and has been shown to reproduce critical obser¬ 
vational constraints such as the low-mass IMF, the brown 
dwarf desert, and the binary statistics of low-mass o bjects 


_ 1 o mei 

dStamatellos fc Whitworthll2009al : lLomax et al.ll2014allm . In 

the third formation scenario mentioned in previous para¬ 
graph, formation by ejection of proto-stellar embryos, ob¬ 
jects that were destined to become Sun-like stars fail to fulfil 
their potential as they are ejected from their natal cloud be¬ 
fore they accrete enough mass to become hydrogen-burning 
stars. 

The presence of discs around substellar objects (and 
associated phenomena, i.e. accretion and outflows) was ini¬ 
tially thought to favour a Sun-like formation mechanism 
(i.e. turbulent fragmentation and collapse of pre-substellar 
cores). However, all three main formation mechanisms pro¬ 
duce substellar objects that are surrounded by discs, albeit 
with different disc fractions. In the turbulent fragmentation 
scen ario substellar object s almost always form with discs 
(e.g. iMachida et al.ll2009ll . Substellar objects that form by 
disc fragmentation also most likely form with discs but these 
discs may be disrupted as these objects are liberated from 
the di sc in which they formed dStamatellos fc Whitworth 
[20093). In tli 6 ejection scenario discs ar e also likely to be 
disrupted but quite a few still survive. iBatd ll2009l') finds 
that at least 10 % of the very-low mass objects formed in his 
simulations have discs with sizes larger > 40 AU. 

Although the presence of discs around substellar objects 
is consistent with all three formation theories, the proper¬ 
ties of these discs may hide clues regardi ng their forma¬ 


tion mechanism. Recently, many authors ([Andrews et al 


2013; iMohantv et al.l l2013l : iRicci et al.l l2014l : [Kraus et al 


2014^ have estimated the masses of discs around young 


low-mass stellar and substellar objects down to a limit of 
~ 10“® ^01 using submillimetre observations . The accre¬ 
tion rates around many low-mass obje cts have also been 
determined down t o 10~^^ ([^atta et all l2004l: 

Calwtet_^lJ_ 2004 1_ Mohantyet_ay 2005^MuxeroUeet^^lJ 


200^ Herczeg&JHffienbra^ 200^ ^itoniucci et al.l 20 111 : 

Rigliaco et al.ll201ll: iBiazzo et aL l2012|j. 


The goal of this paper is to compare these obser¬ 
vations with the theoretical predictions of the disc frag¬ 
mentation model. This is particularly topical as the dis¬ 
covery of many planetary-mass objects at wide separa¬ 
tions (a few tens to a few hundred AU ) from their 


Marois et al.ll2008l; iFahertv et al.ll2009l: 

Ireland et all 

20 111; 

Kuzuhara et al. l201ll. 20131; Aller et al. 

20131; Bailev et all 


2013 

: iRameau et al.ll2013l; iNaud et al.ll2014l; iGalicher et all 

2014 

) has renewed the debate whether these objects have 


formed by core accretion or by fragmentation in the discs 
of their parent stars, or have for med otherwise and were 
later captured by the parent stars dPerets fc KouwenhovenI 
l2012|j . It is also uncertain whether such companions may 
have formed differently than field objects. More wide-orbit 
substellar objects are bound to be discovered with focused 
surveys looking for giant planets (Ge mini Planet Im ager, 
iMacintosh et al.) 20 14; SPHERE/ VLT, iBeuzit et al.l 2008; 
HiCIAO/SUBARU. ISuzuki et all 20091 and therefore their 
properties and the properties of their probable discs may 
be better determined in the near future, providing tighter 
constraints for theoretical models. 

In this paper we present the predictions of the disc 
fragmentation model regarding the masses of discs around 
low-mass stellar and substellar objects (brown dwarfs and 
planets) that are either companion to higher mass stars or 
free-floating. We also determine the accretion rates onto low- 
mass objects and compare them with observations. In Sec¬ 
tion [5] we briefly review the hydrodynamic simulations that 
we use for this study, and in Section |3] we discuss how we 
compute the evolution of the discs around brown dwarfs 
and planets, after these discs have separated from the discs 
of their parent stars. In Section U we present the results of 
the model regarding the disc masses of low-mass objects and 
discuss how they fit with observations, and in Section [S] we 
discuss the accretion rates onto low-mass objects. Finally in 
Section [3 we summerize the main results of this work. 


2 SIMULATIONS OF THE FORMATION OF 
WIDE-ORBIT PLANETS AND BROWN 
DWARFS BY DISC FRAGMENTATION 


2.1 Overview 


The properties of the low-mass stellar and sub-stellar ob¬ 
jects (planets, brown dwarfs, and low-mass hydrogen burn¬ 
ing stars) formed by disc fragmentation have been stud¬ 
ied in detail by Stamatellos e t al. in a series of pa¬ 
pers ([ stamatel los et al.l l2007al ; IStamatellos fc Whitworth! 
l2009al la 2011^. In this paper we use the resu lts of the sim¬ 
ulations of Stamatellos fc WhitworthI ll2009ah to determine 
the properties of the discs around wide-orbit planets, brown- 
dwarfs and low-mass stars that form in the discs of Sun-like 
stars, and the accretion rates onto these objects. 


2.2 Initial Conditions 

IStamatcllos fc WhitworthI ll2009ah performed 12 simulations 
of gravitationally unstable discs around Sun-like stars. These 
simulations are different realisations of the same star-disc 
system, i.e. the properties of the system are the same in all 
simulations; the only difference is the random seed used to 
construct each disc. The star has an initial mass of M* = 
0.7 Mq. The disc around it has an initial mass of Md = 
0.7 Mq and a radius of Ro = 400 AU. The surface density 
of the disc is 


E„(R) 


0.014 Mq / R 

AU^ \A\j) 


( 1 ) 
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and its temperature 


2.4 Results 


r„(R) = 250K ( A 


- 1/2 


+ 10 K. 


( 2 ) 


The disc has an initial Toomre parameter Q ~ 0.9 and there¬ 
fore it is gravitationally unstable by construction. In a re¬ 
alistic situation the disc forms around a young protostar 
and grows in mass by accreting in f alling material from the 
envelope fe.g lAtt wood et al.l 1 20091 : IStaniatellos et allI 2 OIII . 
|2012| : iLomax et al. 2014ah . The disc fragments once it has 


grown enough to become gravitationally unstable at distance 
~ 100 AU from its parent star and this happen before it 
can rea ch the mass assumed by l^amatellos fc WhitworthI 
ll2009ah . In fact even discs with masses ^ 0.25 and 
radii 100 AU can fragment (IStamatellos et ^1201 ill . Such 
disc masses are comparable to the o bserved disc masses in 
young (Class 0, Class I) objects (e.g. lfcrgensen et"^ 20091 : 


iTobin et aLll20l3 : iMurillo et aL|[2013 : Favre et aLll2014h . In 


any case, any evolutionary period with such a massive disc is 
short-lived as the disc quickly (within a few thousand years) 
fragments. 

The large disc m a ss an d size assumed by 
IStamatellos &: WhitworthI (l2009al l ensure that more 
low-mass objects form in the disc to improve the statis¬ 
tical analysis of the results, but the properties of these 
objects (mass, disc mass, disc size) are similar to th e 
ones formed in lower mass discs (IStamatellos et ^l201l|j . 
This is because the characteristic initial mass of objects 
formed by disc fragmentaton is set b y the opacity limit, 


19761: lR,eesl 19761: SilkI 19771: BossI 1 19881: iBovd & Whitworth 

2 OO 5 I: IWhitworth & Stamatellos! l2006l: 

Bolev et al.l I 2 OIOI: 

Kratter et al.ll2010l:lForgan & Rice 201ll: 

Rogers & Wadslev 


formed by fragmentation is the same for lower and higher 
parent disc masses. The parent disc mass (in lower mass 
discs) is distributed among fewer objects and therefore the 
masses of these objects and the masses of their discs are 
similar to the ones that form in higher mass discs. 

The simulations that we use start off with already 
formed discs; therefore disc loading and other interac¬ 
tions with the star forming cloud (which may lead to 
non-axisymmetric discs) are ignored. Simulations that 
take these effects into account (e.g. iTsukamoto et aP 
2^1!^^ have given similar result s to the simulations of 


Stamatellos fc WhitworthI (l2009al l used in the present pa¬ 


per. We therefore do not anticipate the choice of the specific 
set of disc simulations to significantly alter the main conclu¬ 
sions of this paper. 


2.3 Numerical Method 

The evolution and fragmentation of the disc of the par¬ 
ent star is followed using the SPH code dragon which 
treats the radiation tra nsport within the disc with the dif- 
fusion approximati on of IStamatellos et al.l ll2007bf )fsee also 
iForgan et al.ll20(^ . The radiation feedback from the par¬ 
ent star is also taken into account. The code uses time- 
dependent viscosity with parameters a = 0.1, /I = 2 q: (Mor¬ 
ris & Monaghan 1997) and a Balsara switch (Balsara 1995). 


The parent disc is unstable and therefore within a few kyr 
it fragments into 5-11 secondary objects. In the 12 simula¬ 
tions a total of 96 objects are formed. Some of them escape 
and others remain bound to the parent sta r at wide orbits 
(see IStamatellos fc WhitworthI l2009al . l201l|j . Most of these 
objects are brown dwarfs (67%; 13 Mj < M < 80 Mj) 
and the rest are low-mass hydrogen burning stars (30%; 
M > 80 Mj), and planets (3%; M < 13 Mj). These mass 
ranges are set by the hydrogen-burning limit (~ 80 Mj), and 
the deuterium-burning limit (~ 13 Mj). Stars can sustain 
hydrogen burning, whereas brown dwarfs can sustain only 
deuterium burning. Planets cannot sustain deuterium burn¬ 
ing. However, there is no reason for gas fragmentation to 
stop either at the hydrogen-burning limit or the deuterium¬ 
burning limit: the minimum mass of an object that forms 
by gas fragmentation is given from the opacity limit for 
fragmentation (~ 1 — 5 Mj). On the other hand, planets 
that form by core accretion may have masses > 13 Mj (e.g. 
iMolliere fc Mordasini|[201^ . In this paper, we use the term 
planet to refer to objects with mass < 13 Mj regardless on 
their formation mechanism. 


About 70% of the secondary objects that form in the 
parent disc are attended by their own individual discs. These 
discs have masse s up to a few tens of Mj and radi i of a few 
tens of AU (^see IStamatellos fc WhitworthI [2009al l . Out of 
these secondary objects with discs we select 34 single objects 
(i.e. they are not in a binary system with another secondary 
object formed in the parent disc but they may still be bound 
to the parent star) for which the properties of the discs can 
be determined (i.e. the discs are nearly Keplerian). The rest 
of the objects either were binaries, or were attended by disc¬ 
like structures whose properties could not be obtained (e.g. 
discs that were perturbed). Almost all of the objects in the 
sample (33 out of 34) are still bound to the parent star al¬ 
beit in most cases at very wide-orbits (see Fig.[T}. Eventually 
many of these will be liberated an d will become field objects 
dStamatellos fc Whitworthll2009c] l. Therefore, in Sections |4] 
&[5l the properties of the discs of these objects and the ac¬ 
cretion rates onto them will be compared with the observed 
properties of objects that are either wide-orbit companions 
to other stars, or field objects. 

Fig. [T] presents the relation between the masses and 
the semi-major axes of the orbits of these objects. Most of 
these objects are brown dwarfs (with a few of them near 
the brown dwarf-planet boundary of 13 Mj) and a few of 
them planets and low-mass hydrogen-burning stars. Low- 
mass hydrogen burning stars tend to be closer to their 
parent stars than brown dwarfs and planets (the brown 
dwar f desert: iM arcv fc Butleil I 2 OOOI: iGrether fc Lineweaverl 


l2006l : ISahlmann et al. 201ll : iMa fc Gel ~ 2014ll . There are 
many brown dwarf companions to Su n-like stars but at 
these tend to be a t wide separations | Kraus_et_aJ 120081; 
^liert^_et_aL 2009|. 2010|: Kraus et al. 2011 : Evans et alJ 


2 OI 2 I : iReggiani fc Mevejl2013l : iDuchene fc Krausll2013l~l . As 

the above types of objects all form by the same mechanism 
in the disc irrespective of their mass we will analyse their 
disc properties collectively. 

Eig. [5] presents the mass of each disc versus the semi¬ 
major axis of its host object. There is no significant corre¬ 
lation between the two. The disc masses are determined by 
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Figure 1. The masses of objects formed by disc fragmentation 
plotted against the semi-major axes of their orbits around the 
parent star. Most of these objects are brown dwarfs (with a few 
of them near the brown dwarf-planet boundary of 13 Mj) and 
a few of them are planets and low-mass hydrogen-burning stars. 
Low-mass hydrogen burning stars tend to be closer to the parent 
star than brown dwarfs and planets. 
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Figure 2. The disc masses around objects formed by disc frag¬ 
mentation against the semi-major axes of their orbits around the 
parent star. There is no significant correlation between the two. 
The disc masses are probably determined by how these objects 
move in the disc of the parent star and accrete mass from it, 
rather than where they form in the parent disc. 


how these objects move within the disc of the parent star 
and accrete mass from it, rather than where they form in 
the parent disc. 


3 THE EVOLUTION OF THE DISCS AROUND 
BROWN DWARFS AND PLANETS 

The hydrodynamic simulations provide the properties of 
the discs around wide-orbit companions to Sun-like stars 


at the time when 70-80% of the parent disc around the par¬ 
ent Sun-like star has been accreted, either onto the parent 
star or onto the low-mass objects that form in the par¬ 
ent disc. This typically happens within 10 — 20 kyr from 
the start of each simulation. By this point the mass of the 
parent disc has been reduced to < 0.01 Mq. Considering 
that the secondary objects that formed in the parent disc 
are on wide-orbits around the parent star, we do not ex¬ 
pect interactions between the parent disc and the secondary 
discs to be important. Additionally, in a cluster environ¬ 
ment they are likely to be disrupted by stellar flybys and be- 


come free-floating obiects(lHeggielll975 

: iKrouoa et al.ll2003l: 

Parker et al. 20091: Parker & Goodwin 

2009 

: Sourzem et al.l 

20091: Malmberg et al.l201ll: Hao et al. 

2013 

). Therefore, we 


assume that at this point (i) that the secondary discs (i.e. the 
discs around the low-mass objects that form in the parent 
disc) have separated from their parent disc, (ii) that they 
evolve independently (i.e. there are no dynamical interac¬ 
tions between them and the parent disc, or other objects 
the form in the parent disc), and (iii) that no further mass 
from the parent disc is accreted onto them. These assump¬ 
tions are not critical as the accretion of additional mate¬ 
rial onto the secondary disc reinforces our conclusions. To 
compare the properties of the discs around these low-mass 
companions with the observed disc properties of companions 
in nearby young stellar clusters (age ~ 1 — 15 Myr) these 
properties need to be evolved in time. As this is not possible 
to be done by hydrodynamic simulations due to the large 
computational cost, we have employed an analytic model of 
viscous disc evolution. 

We ignore any disc clearing due to photo-evaporation 
fro m radiation from the l ow-mass object hosting the disc 
(see I Alexander et al.ll2oi^ . and references therein). Photo¬ 
evaporation of discs around low-mass obje cts (S: 0.15 Mq) 
could happen (e.g. [Alexander et al.l l2006l l but because of 
our limited knowledge on how UV and X-ray emission from 
low-mass objects would affect their discs, it is difficult to 
ascertain how important photo evaporation is for disc dis¬ 
persion. 

The analytic model we employ assumes the disc (around 
a secondary object) is geometrically thin and evolves vis- 
cous ly under the influence of th e central object’s gravity 
(e.g. iLvnden-Bell &: Pringlel[l974| j. which in this case is the 
planet or brown dwarf (represented as point masses in the 
model). The surface density of such a disc E(i?, t) at polar 
radius R and time t, evolves as follows 


dt 


1A 

RdR 




( 3 ) 


where v{R,t) is the kinematic viscosity l|Prin ^Il98lh . In 
this equation (and in subsequent equations) t — 0 corre¬ 
sponds to the time where these discs are decoupled from 
their parents discs (i.e. the end of the hydrodynamic sim¬ 
ulations). Assuming that the viscosity is independent of 
time and can be expressed as a power law in R, oc R^', 
then the above evolution equation has a similarity s olution 
llLvnden-Bell fc Pringle|[l97^ : iHartmann et al.lll998ll 


E(R,t) 


Md(0)(2 - 7) _(5/2-7)/(2-7) 
27V Rlr-f 



( 4 ) 


where r = R/Ro (Ro is the radius within which 60% of the 
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disc mass is contained initially), and 

T = (5) 

where 

U = Rl/[i{2--tfu{Ro)]. (6) 

The accretion rate onto the central object is then 

• ^ Md(0) (5/2-7)/(2-7) 

* 2{2-y)tJ 

and the disc mass 

Md(t) = . 

It has been argued that obser vations of the discs o f 
T Tauri stars suggest that 7 ~ 1 llHartmann et al]ll998l l 
(i.e. V oc R), and therefore we will adopt this value in the 
present study. The choice of 7 is not critical for the conclu¬ 
sions of this paper. 

We use the Q-viscosity parametarisation 
dShakura fc Sunvaevlll97l3 l 

v = aCsH , (9) 


( 7 ) 

( 8 ) 


where Cs is the sound speed in the disc, H is the disc scale- 
height, and a the viscosity parameter. Assuming that the 
disc is locally vertically isothermal we obtain H = Cs/fl(i?), 
which when used in Eq. and assuming Keplerian rota¬ 
tion, i.e. il{R) = {GMt,/R? f^'^, gives 


ucxaTd . 


( 10 ) 


Using Eq. (fTUl) in Eq. ([^1) and assuming Td{R) oc R 
(consistent with 7 = 1 ) we obtain 


U = 8x10" 


10 - 


-D Ro 

VlOAU 


M* 


1/2 


Y Td 


O.SMgj VlOK 


m 


where Td is the disc temperature at 100 AU. 

We can therefore calculate the disc mass and the accre¬ 
tion rate onto the central object that hosts the disc (planet 
or brown dwarf) at any given time, using the initial disc 
mass Md(0), obtained by the SPH simulations, and using 
Eqs. (O,©, and (na. 


4 THE MASSES OF DISCS AROUND 

LOW-MASS STELLAR AND SUBSTELLAR 
OBJECTS 

Observations of dis c masses (e.g [Andrews et al.l l2013l : 
iMohantv et al.l[2013f l over a wide range of host stellar and 
substellar masses from intermediate mass stars to planetary- 
mass objects suggest a linear correlation between object 
mas s and disc mass , i.e. M a oc M*. 

[Andrews et al.l ([2013[ ~) using 3 different evolutionary 
models for calculating stellar masses they find that stel- 
laiQ mass scales almost linearly with the disc mass, 
Md « 10'"where k = -2.3 ± 0.3, -2.7 ± 0.2, -2.5 ± 
0.2, an d A = 1.4 ± 0.5,1.0 ± 0.4,1 .1 ± 0.4, when us¬ 
ing the iD’Antona &: Mazzitellil ([l997fl (hereafter DM97), 


^ In this context the terms star and stellar are used to refer to 
any objects formed by gravitational instability, therefore includ- 
ing brown dwarfs and p lanets, as well as hydrogen-burning stars 
llWhitworth et ©[200711 . 


[Baraffe et al.l l[l998ll (BCAH 98) and Siesset_aL ([2000[l 
(SDFOO) models, respectively. [Mohantv et al.l ( 20131 1 follow 
a similar approach using the SDFOO models for stars with 
mass > 1.4 Mq, the BCAH98 model for st ars with masses 
0.08 — 1.4 Mq and the dusty models of IChabrier et al.l 
|[2000ll for stellar masses < 0.08 Mq, and si milarly find that 
Md ~ 10”^'^ M*. We note however that both [Andrews et al.l 
|[2013ll and [Mohantv et al.l l[2013ll have assumed that the 
scatter in disc mass is constant for a ll objects irrespective of 
their mass; this may not be the case ([Alexander fc Armitaed 

[ 2 OO 6 II . 


It is evident ( see Fig. 9 in lAndrews et all 2013 and Fig. 9 
in [Mohantv et all 2013) that (i) there is a considerable scat¬ 
ter in the Mdisc — M, relation, and (ii) there are only a 
few definite detections of discs around sta rs with masses 
< 0.1 Mq . For example in the sample of [Andrews et al.l 
|[2013ll using the DCAH98 model there are just 15 definite 
disc detections around stars with mass <0.1 Mq (for 42 
objects only upper limits for the disc masses were derived; 
these upper limits vary from 6 x 10“"* to 1.3 x 10“^ Mq). 
The large scatter in the data points and the small num¬ 
ber of data points at low masses cast doubt to the sugges¬ 
tion that there is a simple relation between stellar and disc 
mass that holds from intermediate-mass stars all the way to 
brown dwarfs and planetary-mass objects. In other words, 
these data do not exclude different scaling relations for low- 
and high-mass objects. Another complication comes from 
the fact that when calculating the relation between stellar 
and disc masses the ages of these objects are not taken into 
account: disc masses are getting smaller with time either due 
to viscous evolution or due to photoevaporation from the 
host star. Thus, considering that discs around low-mass ob¬ 
jects have masses that are low and near the detection limits 
of current observational facilities, it is more likely to observe 
them when they are still young (and therefore have more 
mass). Therefore, it may be expected that the discs around 
brown dwarfs and planets are more massive than what a sim¬ 
ple extrapolation from the Mdisc — M* relation for higher 
mass stars would suggest. The exact effect that the ob- 
ject ages have on the analyses of [Andrews et all l[2013h and 
I Mohantv et al.l ([2013f l is difficult to estim ate as stellar ages 
cann ot be determined accurately enough dSoderblom et all 

[2OI3II . 


The disc masses of the objects formed by d isc fragmen¬ 
tation in the IStamatellos fc Whitwor^ ([2009a[ l simulations 
are plotted against the masses of the objects in Fig. [31 The 
disc masses are calculate d from the disc masses in th e hydro - 
dynamic simulations of IStamatellos fc Whitworth! ([2009ah 
assuming that the discs evolve viscously and using Eq. (©I 
with a = 0.01. The time for which each disc is evolved is cho¬ 
sen randomly between 1 — 10 Myr, so as to emulate the age 
spread of o bserved discs. T h e sam e figure shows the relations 
derived by [Andrews et all ([20131 1 using different evolution¬ 
ary models to calculate the masses of the host star (DM97, 
BCAH98, and SDFOO as marked on graph). 

There is scatter in the calculated disc masses of objects 
formed by disc fragmentation due to differences in the initial 
disc masses (i.e. the mass they have when they separate 
from the disc of the parent star) and their ages. Most of 
these discs are more massive than expected from the scaling 
Mdisc — M* relation (which is mainly determined by higher 
mass stars) by more than a order of magnitude in a few cases. 
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Figure 3. Disc masses of objects formed by disc fragmentation 
versus the masses of the host objects (black squares). The disc 
masses are calculated from the disc masses in the hydrodynamic 
simulations of Stamatellos et al. (2009a) assuming viscous disc 
evolution and a = 0.01. The time for which each disc is evolved 
is chosen randomly between 1 — 10 Myr, so as to emulate the age 
spread of observed discs. On the graph we also plot the best fit 
line (solid black line) and the ±1 (t region from the best fit (dot¬ 
ted black line). The three coloured lines correspond the scaling 
relations derived by Andrews et al. 2013 using different evolu¬ 
tionary models (as marked on graph). The difference between the 
two relations indicate a different formation mechanism for low- 
mass objects. The difference is more pronounced at the extreme 
low-mass regime. 


Additionally, there is no significant dependence between disc 
mass and stellar mass in contrast with higher-mass systems; 
we find a relation log(Md) = —3.7 — 0.005 log(M*) with a 
standard deviation of ct = 0.27. 

Both of the above characteristics are consequences of 
formation by disc fragmentation. When a low-mass object 
forms from gas condensing out in the parent disc, its proper¬ 
ties (and its disc properties) are initially similar to an object 
that forms from a collapsing core in isolation. However as 
this object/disc system moves within the parent disc (but 
before it separates from the parent disc) it accretes more 
gas, and therefore its mass increases. This mass is initially 
accreted onto the object’s disc and then slowly flows onto 
object. Therefore, when a young object that has formed by 
disc fragmentation separates from its parent disc and evolves 
independently, has a more massive disc than it would have 
if it had formed in isolation in a collapsing core. 

This scenario is con sistent with the obse rvations of 
[Andrews et al.l ll2013h and iMohantv et al.l ll2013l) : at least a 
few discs around young low-mass objects are more massive 
than expected. In their samples the detection limit is around 
~ 10“® Mq and a few of the low-mass objects that they ob¬ 
served either have lower-mass discs or no discs at all. These 
may be objects that have either formed by the collapse of a 
low-mass pre-(sub)stellar core like Sun-like stars, or objects 
that have formed by disc fragmentation but have lost their 
discs (through evolution with time, see Fig. [4] or due to in¬ 
teractions within the disc). The presence of low-mass discs 


Figure 4. The evolution of the disc masses of the low-mass ob¬ 
jects formed by disc fragmentation in the simulations of Stamatel¬ 
los et al. (2009a). The disc masses are calculated assuming viscous 
disc evolution with a viscosity parameter a = 0.01. Disc mass de¬ 
creases with time. Due to the difference in the initial disc masses 
there is a wide range of disc masses for each given time. 


around low-mass objects are consistent with both formation 
scenarios but the presence of relatively high-mass discs are 
indicative of formation by disc fragmentation. Observations 
of disc masses around very low-mass objects (< 10 Mj), 
where the predicted Mdisc — A7* relation for young objects 
diverges significantly from the established Mdisc — M« re¬ 
lation derived for higher mass stars, will further test the 
disc fragmentation model. ALMA has the required sensitiv¬ 
ity and spatial reso l ution to observe such small discs. For ex¬ 
ample, [RiccOri^ ll2ni4h . have estimated disc masses down 
to ~ 0.8 — 2.1 Mj in three young low-mass objects in the 
Taurus star forming region. 


5 ACCRETION RATES ONTO WIDE-ORBIT 
LOW-MASS OBJECTS 


The accretion rates onto low- and higher-mass objects may 
also relate to their formation mechanism. In some cases 
it is possible to derive accretion rates even when the disc 
that provides the material for accretion is not detectable 
in the sub-mm, where disc masses are usually measure d 
llHerczeg et al.ll2 009l:ljoergens et ^l2013l:IZ hou et ^|2014|) . 

For example, Herczeg et al. ( 20091 ) and IZhou et al.l ( 20141) 


estimate the accretion luminosity from the excess line and 


continuum emission; for low-mass objects 
accretion rates down to ~ 10“^® Mq yr“ 
It has been argued that, similarly to ■ 
lation mentioned in the previous section, 
between accretion rate onto a star and its 
suggested that this relation holds from 
stars down to brown dwarfs, namely that 
a ^ 1.0 — 2.8, albeit with a large scatter ( 

they can estimate 

he Mdisc — M* re- 
there is a relation 
mass. It has been 
ntermediate-mass 
M* oc M“, where 

Natta et al. 2004: 

Calvet et ahl 20041: Mohantv et al.l 2005 

: iMuzerolle et al.1 

20051: iHerczeg & Hillenbrandl 20081: lAntoniucci et ahl 201 ll: 

Biazzo et al. 20121). 
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Figure 5. The accretion rates onto stars against their masses 
for a wide range of stellar masses. The black squares correspond 
to the objects formed by disc fragmentation in the simulations of 
Stamatellos et al. (2009) with accretion rates calculated assuming 
viscous disc evolution with a = 0.01. The time for which each disc 
is evolved is chosen randomly between 1—10 Myr, so as to emulate 
the age spread of observed discs. On the graph we also plot the 
best fit line (solid black line) and the dilcr region from the best fit 
(dotted black line). The remaining points (coloured) correspond 
to observational data as marked on the graph. On the graph we 
also plot the best fit line for the observations (solid blue line) and 
the dilcr region from the best fit (dashed blue lines) as estimated 
by Zhou et al. (2014). 


The accretion rates onto stars for a wide range of stel¬ 
lar masses are plotted against the stellar masses in Fig. O 
The accretion rates shown here have all been measured di¬ 
rectly from excess Ba l mer continuum emission in the U- 


band (iGullbrinEf et alJ 19981: Herczeff Kr. Hillenbrandl 20081: 

Herczeff et al.ll20od:lRle'liaco et al.ll201ll.l2012l: Inelebv et al.l 

20131: Alcala et al.ll20l4 IZhou et all 2014 

). In the same fig- 

ure the best fit line that was calculated by 

Zhou et al.l (12014) 


is also plotted. It is evident from the graph that there is con¬ 
siderable scatter in M* — M* relatio n, that may reflect a dif¬ 
feren c e in the disc initial cond itions ijAlexander Armita^ 
l2006l : IPullemond et al ][200i). A part of the scatter could 
also be attributed to the different ages of the systems plot¬ 
ted in Fig. m accretion rates drop as stars age (see Fig. [6]). 

The detection limits of accretion rates are relatively low 
for planets and brown dwarfs. Most objects with excess emis¬ 
sion in the IR also have measured U-band accretion rates; 
thus it is expected that there is no bias towards detect¬ 
ing only younger objects with higher accretion rates. In fact 
most of the observed objects exhibit low accretion rates. The 
estimated accretion rates for most of the low-mass objects 
(< 0.1 Mq) are consistent with the M* — M* scaling rela¬ 
tion derived for higher-mass stars. In fact in a few cases the 
accretion rates are lower than expected. However in a few 
cases, like the thr ee planetary-mass companions observed by 
IZhou et al.l (120141 ') (GSC 06214-00210 b, GQ Lup b, and DH 
Tan b) the accretion rates are higher than expected; these 
accretion rates are an order of magnitude higher than what 
is expected from the M* — M* relation. 


10^=5 
10-6 
10 -^ 
y 10-8 
0 10-6 

10 -" 

10-12 
10-18 

0 2 4 6 8 10 

t (Myr) 

Figure 6. The evolution of the accretion rates of the objects 
formed by disc fragmentation in the simulations of Stamatellos et 
al. (2009a) . These accretion rates are calculated using the viscous 
evolution model (Eq. l[7]l with a = 0.01). There is a wide range of 
accretion rates for a specific age due to the spread in the initial 
disc masses. The three red triangles corresponds to the obser¬ 
vations of Zhou et al. (2014). Considering the large uncertainties 
(~ 1 —5 Myr) in the estimated ages these relatively high accretion 
rates are consistent with the predictions of the disc fragmentation 
model. 


In Fig. [5] we also plot the accretion rates of the ob- 
jects formed by disc fragmentat ion in the simulations of 
IStamatellos fc Whitworth[l|2nn9ah . These accretion rates are 
calculated using the viscous evolution model (Eq. ([7}) with 
a = 0.01. The time for which each disc is evolved is cho¬ 
sen randomly between 1 — 10 Myr, so as to emulate the 
age spread of observed discs. There is no significant correla¬ 
tion between the accretion rate and the mass of the object; 
we find a relation log(M,) = —10.5 — 0.12 log(M*), with a 
standard deviation of cr = 0.3. Moreover, in a few cases the 
accretion rates are higher than expected from the M* — M« 
scaling relation. In the model that we present here, this is 
due to the higher initial mass of the discs of these objects. 
As mentioned in the previous section, these secondary discs 
grow in mass as they move within the discs of their parent 
stars (before they start evolving independently). Therefore, 
we suggest that the relatively high accretion rates are indica¬ 
tive of formation by disc fragmentation. On the other hand, 
low accretion rates are consistent with both formation by 
disc fragmentation or formation by the collapse of low-mass 
pre-(sub)stellar cores. In the former case low accretion rates 
could be due to time evolution (accretion rate drops with 
time; see Fig. El or due to disruption by interactions with 
other objects in the parent disc. 

Observations of acc retion rates aroun d very low-mass 
objects (;$ 10 Mj; e.g. IZhou etTll2014h . where the pre¬ 

dicted M* — M* relation relation diverges significantly from 
the established M, — M* relation relation derived for higher 
mass stars, will further test the model presented here. 
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6 THE EFFECT OF THE VISCOSITY OF 

SECONDARY DISCS 

We have so far assumed in our analysis that the physical pro¬ 
cesses for redistributing angular momentum are the same 
for discs of T Tauri stars and for discs of lower mass ob¬ 
jects (brown dwarfs, planets). However, this may not be the 
case. It has been argued that the magneto-rotational in- 
stabil ity may not be effective in discs around low-mass ob¬ 

jects llKeith fc Wardlell20l4 ISzulagvi et al.l[20l3 : iFuiii et ^ 
|20141 . which means that the effective viscosity in such discs 
should be smaller than the one presnmed for T Tanri star 
discs (q = 0.01). However, these studies have focused on 
discs around Jovian planets on Jovian orbits, i.e. orbits re l- 
atively close to the central stars (e.g. iGressel et al.l[201^ . 
In onr study we focus on wide-orbit low-mass companions 
(see Fig. [21, whose discs are more extended as they not lim- 
ited by the Hill radii of their h ost secondary objects (see 
IStamatellos fc Whitworthll2009a| l. These discs could be mas¬ 
sive enough so that angular momentum can be effectively 
transported by gravitational torques. 

Nevertheless, our knowledge of the effective viscosity in 
such discs is limited, and it is important to examine the ef¬ 
fect that the assumed disc viscosity has on the conclusions 
of onr study. In Figs. [7] and [8] we present the predictions 
of our model for low-viscosity discs {a = 0.001) and for 
high-viscosity discs (a = 0.05). As expected, low-viscosity 
discs evolve slower and their masses and accretion rates re¬ 
main higher for longer. Therefore, in this case the differences 
between the predicted M^isc — M, and M* — M* relations 
for disc fragmentation and the observed relations for higher 
mass stars are more pronounced (see black lines in Figs. [T] 
d. The opposite holds for high-viscosity discs {a = 0.05; see 
brown lines in Figs. izld. 


7 CONCLUSIONS 

We suggest that substellar (planetary-mass objects and 
brown dwarfs) and low-mass stellar objects (low-mass hy¬ 
drogen burning stars) that form by disc fragmentation, have 
disc masses and accretion rates that (i) are independent of 
the mass of the host object, and (ii) are higher than what is 
expected from scaling relations derived from their interme¬ 
diate and higher-mass counterparts. These low-mass objects 
form similarly to higher-mass objects by self-gravitating gas 
but as they move within the gas-rich parent disc their in¬ 
dividual discs accrete additional material; therefore before 
these objects separate from their parent discs and evolve in¬ 
dependently (i.e. within a few kyr), their discs grow more 
massive and the accretion rates onto them are higher than 
if they were formed in isolation in collapsing low-mass pre- 
(sub)stellar cores. The assumption of independent evolution 
is not critical as if these secondary discs were still interacting 
with their parent disc they would accrete additional mate¬ 
rial reinforcing the above conclusion. However, we do not 
expect additional accretion to be important. 

Observations of disc masses and accretion rates of low- 
mass objects are consistent with the predictions of the disc 
fragmentation model. Although the presence of low-mass 
discs (or lack of discs) and low accretion rates (or no accre¬ 
tion at all) may be attributed to disc evolution and/or disc 



Figure 7. Disc masses of objects formed by disc fragmentation 
versus the masses of the host objects for different disc viscosities. 
The best fit lines are calculated similarly to the ones in Fig. (4] 
assuming that a = 0.001 (black solid line) or a = 0.05 (brown 
solid line). The dotted lines correspond to the iltr region from the 
best fit for each case. The other three coloured lines correspond to 
the scaling relations derived by Andrews et al. 2013 (see Fig. 11. 
The differences between this relation and the ones derived in this 
paper (i.e. for objects formed by disc fragmentation) are more 
pronounced for low-viscosity secondary discs. 



Figure 8. The accretion rates onto stars against their masses for 
a wide range of stellar masses. The best fit lines are calculated 
similarly to the ones in Fig. (5] assuming that a = 0.001 (black 
solid line) or a = 0.05 (brown solid line). The dotted lines cor¬ 
respond to the dilcr region from the best fit for each case. On 
the graph we also plot the best fit line for the observations (solid 
blue line) and the ilcr region from the best fit (dashed blue lines) 
(Zhou et al. 2014). The differences between this relation and the 
ones derived in this paper (i.e. for objects formed by disc fragmen¬ 
tation) are more pronounced for low-viscosity secondary discs. 
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disruption due to interactions with other objects within the 
parent disc, relatively high disc masses and high accretion 
rates are suggestive of formation due to disc fragmentation. 
We therefore suggest that low-mass objects that have discs 
with masses higher than expected (or equivalently accretion 
rates onto them higher than expecte d), such as GSC 0 6214- 
00210 b, GQ Lup b, and DH Tau b (IZhou et aLll20l4) . are 
young objects that have formed by disc fragmentation. 

The disc fragmentation model can further be tested 
by observations of disc masses and accretion rates of very 
low-mass objects (J: lOMj). At these very low-masses the 
Mdisc — and M* — M* relations predicted by the model 
presented here diverge significantly from the corresponding 
relations established for higher-mass stars. We suggest that 
future analyses of the Mdisc — M« and M* — M« relations 
should separate the sample into two subgroups, low-mass 
(< 0.2 Mq) and higher-mass (> 0.2 Mq) objects, so as to 
test whether these objects obey different scaling relations. 

The intense interest in wide-orbit and free-floating plan¬ 
ets has given momentum to the development of instruments 
with high sensitivity and good spacial resolution. Therefore 
observations in the near future are expected to deliver many 
more such low-m ass objects. ALMA is already deliv ering 
such observations (iRicci et al.ll20lJ : iKraus et al.ll20l3) . The 
study of these objects, their disc properties and the accre¬ 
tion rates onto them (if they are still young) will provide 
further constraints regarding their formation mechanism. 
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